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Abstract 


Objective: To investigate the effect of long-term Tai Chi Chuan (TCC) practice on practitioners’ brain functional 


specialization compare with the TCC novices. 


Methods: 22 TCC practitioners (52.4 + 6.8 years; 7 males; educated years:12.18 + 3.03 years) and 18 healthy controls 


(54.8 + 6.8 years; 8 males; education years:11.78 + 2.90 years) matched by age, sex, and education were enrolled. 
Participants underwent functional magnetic resonance imaging (fMRI) scanning and cognitive test to measure the 
differences in functional specialization and cognitive function. Functional specialization was evaluated by voxel-mirrored 
homotopic connectivity (VMHC) method. 

Results: Lower middle frontal gyrus (MFG) VMHC in TCC practitioners compared to controls. For TCC practitioners, the 
longer they practice, the lower their VMHC in precentral and precuneus. TCC practitioners showed better cognition 
performance. 

Limitations: Relatively small sample size. 

Conclusions: Changed VMHC indicated that TCC practice could enhance functional specialization in the middle frontal 
cortex of practitioners, which may be associated with higher-order cognitive ability. 
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1. Introduction 

Increasing amount of evidence indicates that mind-body exercises (e.g., Yoga, Tai Chi Chuan 
(TCC), and certain Health-Qigong forms) are not only efficient relaxation skills promoting good 
health [1], but can also act as alternative treatments for symptomatic management in individuals with 
a number of medical conditions [2-4]. TCC is a typical form of mind-body exercise originating in 
ancient China, which involves meditation, mental concentration, musculoskeletal stretching, and 
diaphragmatic breathing with choreographed movements [5, 6]. Given increasing evidences 
demonstrate its beneficial effect on behavioral performance, it has been accepted for symptomatic 
management of some chronic diseases. For instance, a recent study by Wang et al. [7] observed that 
12-month TCC intervention resulted in similar or greater improvement in disease specific symptoms 
in patients with fibromyalgia compared to those engaging in the same dose of aerobic exercise. Some 
evidence also indicates that TCC practice can regulate emotion and relieve the symptoms of mood 
disorders [8, 9]. Additionally, a number of behavioral studies demonstrated that TCC can improve 
cognitive function (e.g., executive function, reaction time) [10, 11]. To better understand the neural 
mechanisms underlying such improvements, magnetic resonance imaging (MRI) studies were 
performed on TCC practitioners. These studies reported that experienced TCC practitioners exhibited 
greater cortical thickness [11] and larger brain volumes, relative to healthy controls [12, 13]. 
Functional MRI (fMRI) studies indicated that long-term TCC practice optimizes functional 
connectivity (FC) in brain regions which play key roles in cognitive-motor integration in adults [14- 
17]. These findings indicate that TCC practice may reshape brain function and could be associated 
with cognitive improvements. 

Brain functional specialization emphasizes hemispheric processing, which allows the expanded 
human brain to minimize between-hemisphere connectivity and distribute domain-specific processing 
functions [18]. The improvement in functional specialization between hemispheres may enhance the 
efficiency of cognitive processing. Currently, there is no direct evidence supporting the effect of 
extensive TCC training on functional specialization. The potential role of functional specialization 
between the two hemispheres in the effect of TCC practice on cognition remains largely unknown. 


Voxel-mirrored homotopic connectivity (VMHC), a method which quantifies the resting-state FC 


between each voxel in one hemisphere and its corresponding voxel in the opposite hemisphere, can be 
used to evaluate the functional specialization between the two hemispheres [19]. VMHC was 
commonly used to investigate the functional homotopic connectivity in cortical regions, including 
areas involved in visual, motor, somatosensory, and higher order association functions in healthy 
individuals [19, 20]. Additionally, this analysis was also applied to investigate the aberrant 
interhemispheric FC in patients with neuropsychiatric disorders, such as schizophrenia, major 
depressive disorder, and panic disorder [21-23]. Therefore, VMHC can be an effective approach for 
the evaluation of functional specialization induced by TCC practice. 

In the present study, we employed fMRI to investigate functional homotopy induced by TCC 
practice, and to examine the neurobiological implications of altered FC. We adopted a cross-sectional 
design to investigate the differences in functional specialization (by measuring VMHC) and the 
differences in behavioral performance (by assessing cognitive task performance) between TCC 
practitioners and healthy controls. It was hypothesized that TCC induces changes in VMHC in brain 
regions responsible for higher-order information processing (prefrontal cortex), due to its demand for 
complex movements. Such changes would provide a neural mechanism for the cognitive outcomes 


induced by mind-body exercise from a functional specialization perspective. 


2. Methodology (Design/Approach) 
2.1 Participants 


Twenty-two TCC practitioners (median age, 52.4 + 6.8 years; 7 males) and eighteen healthy 


controls (54.8 + 6.8 years; 8 males) matched by sex, age, and education were enrolled in this study. 
TCC practitioners were recruited from local TCC activity centers. On average, participants had 14 + 8 
years of TCC experience. The practice frequency (number of practice sessions per week) and the 
duration of each session were used to calculate the total hours of TCC practice per week, which served 
as an index of the intensity of TCC practice. TCC styles mainly included Yang, Wu, Sun, and modified 
Chen. The length of formal TCC practice ranged from 30 to 90 minutes per session, with the majority 
of TCC practitioners (85%) attending daily sessions. Healthy controls were recruited from the local 


community and had no previous experience of regular physical activity and mind-body exercise (e.g., 


TCC, Qigong, and Yoga). Each participant received financial compensation for their time. All 
participants met the following criteria: (1) right-handed; (2) not claustrophobic; (3) no metal implants 
in the body; (4) no mental or neurological disorder; and (5) no history of alcohol or substance abuse. 
Written informed consent was obtained from all participants before the experiment. This study was 
approved by the Institutional Review Board and followed the ethical standards described in the 
Declaration of Helsinki throughout the study. Before scanning, all participants completed screening 
forms of the MRI center to exclude any participant ineligible for fMRI. 

2.2 Behavioral Assessment 

Attention Network Test (ANT) was originally developed by Fan et al. [24] to quantify different 
dimensions of attention, including alerting, orienting, and executive attention. This test has been 
widely used to investigate multiple attentional dimensions across different populations with good 
reliability and validity [25]. Due to time limitations, 10 participants from each group (a total of 20 
participants), matched by age, sex, and years of education, completed ANT before undergoing fMRI 
scanning in this study, in order to measure executive attention. 

The presentation of stimuli and the collection of behavioral data were performed using E-prime 
Software 2.0. One target stimulus and four flanker stimuli as a stimulus sequence were presented in a 
horizontal row. The target stimulus was an arrow pointing either to the left or to the right and two 
flanker stimuli were arrows or neutral stripes which flanked the target on each side. The task involved 
three conditions: congruent, incongruent, and neutral. In the congruent condition, all arrow stimuli 
pointed in the same direction. In the incongruent condition, the target arrow pointed in the direction 
opposite to flanker arrows. In the neutral target condition, all four flanker stimuli were stripes. A 
fixation cross was presented in the middle of the screen, with the stimulus sequence presented at a 
visual angle of 1.1 degrees above or below it. Before a stimulus sequence was presented, the location 
and timing were cued in four different ways: (1) an asterisk was presented at the location of the 
fixation cross, providing the location information of target stimulus; (2) two asterisks appeared 500 ms 
before target onset at the fixed location of 1.1 degrees of visual angle above and below the center of 
the screen, providing information about the timing of the target stimulus; (3) an asterisk was presented 


at the future location of the target stimulus 1.1 degrees of visual angle above or below the center of the 


screen, providing both the timing and the location information of stimulus sequence; and (4) no cue 
was given, such that participants could not obtain any timing or location information regarding the 
upcoming stimulus sequence. 

A complete ANT consists of one practice block with 24 trials and three test blocks with 96x3 
trials. The test takes approximately 25 minutes to complete. A short break was provided between the 
blocks. In each trial, a fixation cross was presented in the middle of the screen from 400 to 1600 ms, 
followed by a cue presented for 100 ms and a central fixation for 400 ms. The stimuli sequence was 
then presented for 1700 ms. If a response was given within 1700 ms, the presented time was shorter. 
Finally, a fixation cross was presented again. Participants were instructed to respond as fast and as 
accurately as possible to the target stimulus. When the target arrow pointed to the left, participants 
were asked to left-click with their left thumb. If the arrow pointed to the right, participants were asked 
to right-click with their right thumb. The reaction time and accuracy were recorded in each trial. The 
final score of executive attention equals the average reaction time in incongruent flanker conditions 
minus the average reaction time in congruent flanker conditions. 

In addition to ANT, all participants were administrated the Number Cancellation Test (NCT) 
[26]. This test was developed to measure selective attention and consists of one practice sheet and four 
test sheets (Cronbach a = 0.94, Test-Retest Reliability = 0.90). Each sheet contains 200 single digits 
with symbols below or above them. Forty-four of these digits were targets, while the rest of them were 
non-targets. Targets were digits “9” with two symbols below, above, or on either side, except those 
right behind digit “5.” Participants were asked to cross out targets as quickly as possible. Only 1 min 
was permitted for each sheet. Between two sheets, they were provided a 10 to 20 s break. The NCT 
score was calculated using the following formula: 

NCT score = hits number - false alarms number x 44/156. 

2.3 MRI Data Acquisition 

Brain images were acquired using a 3T Trio Tim Siemens scanner with a 12-channel head matrix 
coil. Participants were instructed to keep their eyes closed, remain relaxed, stationary, and awake 
during the resting scans. Using an echo planar imaging sequence (TR = 2000 ms, TE = 30 ms, FA = 


90°, slice thickness = 3.0 mm, field of view = 200 x 200 mm”, voxel size = 3.4 x 3.4 x 4.0 mm), 243 


volumes with 33 axial slices were obtained from each participant. For better registration and overlay 
of brain activity, high-resolution structural images were also acquired through three-dimensional 
magnetization prepared rapid gradient echo (3D MPRAGE) sequence (TR = 2530 ms, TE = 3.39 ms, 
FA = 7°, slice thickness = 1.33 mm, voxel size = 1.33 x 1.0 x 1.33 mm). 
2.4 Image Preprocessing 

The preprocessing of fMRI data was conducted using a MATLAB toolbox Data Processing 
Assistant for Resting-State fMRI (DPARSF) [27], which is based on functions in Statistical 
Parametric Mapping (SPM) and Resting-State fMRI Data Analysis Toolkit (REST). To ensure the 
steady-state of longitudinal magnetization, the first 10 images of each participant’s scan were 
removed. Slice timing and realignment were then conducted to correct for temporal differences and 
head motion. Next, T1-weighted images were co-registered to the mean functional volumes and 
segmented into gray matter, white matter, and cerebrospinal fluid. Subsequently, the 24 head-motion 
parameters (6 head-motion parameters of current time point, 6 head-motion parameters of previous 
time point, and 12 corresponding squared items), white matter signal, cerebrospinal fluid signal, and 
the global signal were regressed from functional volumes which were then spatially normalized to the 
Montreal Neurological Institute (MNI) space and resampled to 3 mm x 3 mm x 3 mm resolution. 
Furthermore, linear detrending and temporal bandpass filtering (0.01 — 0.08 Hz) were performed to 
reduce the effects of low-frequency drift and high-frequency physiological noise. Finally, functional 
volumes were registered to a study-specific and symmetric MNI template which was automatically 
generated by DPARSF based on data from all study participants and smoothed (FWHM = 4 mm) to 
improve the space symmetry of each voxel between two hemispheres, increasing the validity of 
VMHC calculations. 
2.5 VMHC Analysis 

VMHC analysis was performed using DPARSF. To minimize the effects of differences between 
the two hemispheres, a symmetric brain template of left-right hemisphere was generated. Next, the 
homotopic functional connectivity was computed as the Pearson correlation between the residual time 


series of every pair of symmetrical interhemispheric voxels. Pearson correlation coefficients were then 


transformed into Fisher's z score to improve normality. The final z-maps were referred to as VMHC 
maps and were used for further statistical analyses. 
2.6 Statistical Analysis 

In behavioral studies, chi-squared tests and two-sample t-tests (two-tailed) were used to compare 
the demographic and behavioral differences, respectively, between TCC practitioners and healthy 
controls using SPSS 19.0 software. Group differences in age, education, and behavioral parameters 
between TCC practitioners and healthy controls were tested by two-sample t-tests. Chi-square test was 
used to determine the difference in sex between the two groups. 

On the neuroimaging level, considering age, sex, and education as covariates, a two-sample t-test 
was performed on VMHC maps using SPM8 to detect the differences in homotopic connectivity 
between TCC practitioners and healthy controls. The statistical significance level for group 
differences in VMHC was corrected with family-wise error approach (p < 0.05) to account for 
multiple comparisons. Furthermore, we calculated the Pearson’s correlation coefficients between 
VMHC values and TCC behavioral results in TCC group to investigate the association between 
interhemispheric functional connectivity and TCC-induced higher order cognitive function. 


3. Results 


3.1 Demographic and Behavioral Data 


The demographic characteristics and behavioral data of participants are described in Table 1. 


There were no significant differences in sex (y?= 0.673, p = 0.517), age (t (38) = 1.149, p = 0.258), or 
education level (t (38) = 20.435, p = 0.666) between the TCC and control groups. Relative to the 
control group, TCC participants showed marginally better performance in NCT (t (18) = 1.951, p = 
0.067). However, no significant group differences were detected both in reaction time (t (18) = -1.227, 
p = 0.236), and accuracy of ANT (t (18) = 1.472, p = 0.158). 

3.2 Functional Brain Imaging Results 


VMHC values in the TCC group were significantly lower than those in the control group in the 


middle frontal gyrus (MFG, MNI coordinates x = +44, y = 41, z= 13; p < 0.001, corrected for family- 
wise error; Figure 1). Correlation analyses demonstrated no significant relationship between VMHC 


values in the MFG and TCC experience (r = 0.264; p = 0.274) or cognitive task scores (ANT reaction 


time, r = -0.366, p = 0.419; ANT accuracy, r = -0.143, p = 0.759; NCT scores: r = 0.248, p = 0.592). 
Furthermore, we evaluated the correlation between TCC experience and VMHC values in all brain 


regions, which demonstrated a significant negative correlation between the total training hours of TCC 


practice and VMHC values in the precentral cortex (MNI coordinates x = +31, y = -23, z= 69; r 


-0.983; p = 0.000; Figures 2A and 2B) and in the precuneus (MNI coordinates x = +18, y = -81, z 


42; r = -0.807; p = 0.000; Figures 2C and 2D). 
4. Discussion 

In this study we utilized the cross-sectional designs to examine if TCC exercise can induce a 
change in functional homotopy, and to elucidate the relevant biological basis of altered brain function. 
To the best of our knowledge, it is the first study to utilize VMHC analysis method to explore regional 
specific difference of functional homotopy between TCC experts and controls. We observed decreased 
VMHC values in MFG in TCC group and significant negative relationship between total TCC time 
and VMHC in precuneus and precentral. These findings indicate that the possibility of local plasticity 
in functional specialization of brain as function of extensive training. 

Left-right asymmetry is a key feature of the human brain structure and function [28]. A largest- 
ever analysis of cerebral cortical asymmetry found that regional asymmetries in cortical thickness 
and/or surface area are extensively present in the inferior frontal gyrus, transverse temporal gyrus, 
parahippocampal gyrus, and entorhinal cortex of healthy individuals [29]. Gotts et al. [30] also pointed 
out that two distinct forms of functional specialization were present in the left versus the right cerebral 
hemisphere. The left hemisphere tends to interact more exclusively with itself, while the cortical 
regions of the right hemisphere involved in visuospatial and attentional processing tend to interact in a 
more integrative fashion with both hemispheres. The widespread lateralization across species may 
suggest an evolutionary advantage associated with the specialization of the hemispheres [31]. An 
earlier study developing the VMHC approach found that homotopic regions in the prefrontal cortex 
(particularly those in the ventrolateral and dorsolateral regions responsible for attention, cognition, 
and language) exhibited lower coherence [19]. This finding indicates that functional specialization, 


especially in the higher-order cerebral cortex, is closely associated with human cognition. 


The most interesting finding in our study is the detection of decreased VMHC in MFG. The 
prefrontal lobe is generally considered to be the high-order cognitive and executive center which 
integrates the information from internal sources and the external environment [32, 33]. As a part of the 
prefrontal lobe, MFG plays an important role in complex information processing related to goal- 
directed cognition [34]. Several studies demonstrated that the alteration in functional activity in the 
MFG is associated with improved cognitive performance, which can be induced by extensive TCC 
training. A six-week TCC intervention elicited increased spontaneous activation in the right MFG 
which was associated with improved cognitive performance. Additionally, there was a significant 
correlation between amplitude of low frequency fluctuations (ALFF) changes in the right MFG and 
changes in executive function among participants who were involved in a six-week intervention 
including Tai Chi exercise [35]. The functional difference in MFG between TCC experts and controls 
possibly suggest greater functional plasticity in MFG. Decreased VMHC in MFG suggests that 
functional specialization of MFG may reflect enhanced utilization of the functional cortex, which 
resulted in better cognitive ability among TCC practitioners in this study. 

Decreased VMHC in precuneus was associated with TCC experience. As a component of the 
default-mode network (DMN), precuneus functions as the brain region processing emotional 
information and multiple higher-order cognitive functions [36, 37], such as attention shifting, memory 
retrieval, and spatial coordination [38]. Increasing evidence demonstrates that meditation practice can 
induce structural and functional changes in the precuneus [39, 40]. A recent TCC study demonstrated 
that TCC experts had decreased fractional ALFF in DMN, compared to healthy controls [17], which 
directly supports the perception that the precuneus brain region can be altered by meditative training. 
Importantly, a number of studies confirmed that meditative exercise can change the precuneus 
asymmetry [41, 42]. The decreased VMHC in precuneus provided further evidence that TCC is likely 
to alter the inter-hemispheric functional integration in the precuneus. 

In addition to the precuneus, precentral sulcus is another brain region which exhibits significant 
association with TCC experience. Located in the primary motor cortex, precentral sulcus is 
responsible for observation and execution of motor tasks [43]. A study reported that the lateralization 


of sensorimotor networks (including precentral gyrus) decreases with aging [44], suggesting that aging 


may lead to decreased hemispheric specialization. Fujiyama et al. [45] found that aging adults show 
bilateral recruitment of primary motor cortex while younger adults show a suppression of the 
ipsilateral motor cortex during motor tasks. It is plausible that TCC practice enhances lateral function 
in the sensorimotor cortex to protect against aging, thereby improving the capability of motor-related 
information processing. Luders et al. [46, 47] suggested that long-term meditation practitioners may 
have younger brains. Since TCC shares a number of mental components with meditation (such as 
concentration and relaxation), these findings may help us better understand the neural mechanism 
underlying the beneficial effect of TCC practice on motor ability. 

One limitation of the current study is the relatively small sample size. In particular, the sample 
size in behavioral testing involved only half of the participants. This may partly explain why 
significant correlation between VMHC in MFG and cognitive performance was not observed in the 
cross-sectional study. Therefore, future studies with large sample population are needed to confirm the 
role of functional homotopy in the precuneus and precentral sulcus. 

To sum up, decreased VMHC in the MFG in TCC participants suggesting that extensive TCC 
practice may strengthen functional specialization in the MFG. The findings extend our previous 
knowledge by suggesting functional homotopy as a biomarker of plasticity induced by mind-body 
exercise. Additionally, although the association between decreased VMHC and cognitive performance 
was not observed in this study, higher cognitive performance in TCC group provides insights for 
further investigations aiming to examine the biological relevance of enhanced functional 


specialization with respect to mind-body practice from a cognitive perspective. 
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Figure Legends 
Fig. 1 Differences in VMHC values between TCC experts and control participants in experiment 1. 


Blue areas indicate decreased VMHC values in MFG in TCC experts. The color bar represents the 


significance of group analysis. 


Fig. 2 The correlation between the duration of TCC practice and VMHC values of all voxels on the 
global level. (A) Location of the precentral. (B) The relationship between total training hours in TCC 
experts and their VMHC values in precentral cortex. (C) Location of the precuneus. (D) The 


relationship between total training hours in TCC experts and their VMHC values in precuneus. 
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Tables 


Table 1. Demographic characteristics and behavioral data of participants 


Characteristic TCC Control Statistic P value 
(two-tailed) 
Age (years) 52.36 + 6.88 54.83 + 6.77 t (38) = 1.149 0.258 
Education (years) 12.18 + 3.03 11.78 + 2.90 t (38) = 20.435 0.666 
Sex (male/female) 7/15 8/10 = 0.673 0.517 
NCT 25.09 + 3.53 21.68 + 4.25 t (18) = 1.951 0.067 
ANT (RT', ms) 98.50 + 31.74 123.01 + 54.63 t (18) = -1.227 0.236 
ANT (AC?) 0.9948 + 0.0060 0.9899 + 0.0086 t (18) = 1.472 0.158 


'RT = Reaction Time; ? AC = Accuracy 


